Hramov A.E., Koronovskii A. A., Popov P.V. Incomplete noise-induced synchronization of spatially 
extended systems. Phys. Rev. E. 77, (2008) 036215 



Incomplete noise-induced synchronization of spatially extended 

systems 

Alexander E. Hramov, 1 Alexey A. Koronovskii, 1 and Pavel V. Popov 1 

1 Faculty of Nonlinear Processes, Saratov State University, 
Astrakhanskaya, 83, Saratov, 410012, Russia 
(Dated: June 27, 2008) 

Abstract 

A new type of noise-induced synchronous behavior is described. This phenomenon, called incom- 
plete noise-induced synchronization, arises for one-dimensional Ginzburg-Landau equations driven 
by common noise. The mechanisms resulting in the incomplete noise-induced synchronization in 
the spatially extended systems are revealed analytically. The different model noise are considered. 
A very good agreement between the theoretical results and the numerically calculated data is 
shown. 
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Introduction 



Noise-induced synchonization [1|, [2j, |3| is an ubiquitous phenomenon in nonlinear science. 
It arises as the interplay between determined and random dynamics 4(, with both the 
synchronization and noise influence being recently the subjects of considerable interest of 
scientific community. Indeed, on the one hand the synchronous behavior of nonlinear systems 

nnnn 

has attracted great attention of researchers for a long time [5|, [6|, 13, |8| . On the other hand 
discovering the fact that fluctuations can actually induce some degree of order in a large 
variety of nonlinear systems is one of the most surprising results of the last decades in the 
field of stochastic processes fly , Moreover, both these phenomena are relevant for 



etiological aad other systems described in *™ of nonlinear dynamics 



[see, e.g., 



13|,[14|,[15(; 



Noise-induced synchronization (NIS) means that the random signal influencing two iden- 
tical uncoupled dynamical chaotic systems u(t) and v(t) (starting from the different initial 
conditions u(t ) and v(t ), u(to) 7^ v (^o)) results in their synchronous (i.e., identical) be- 
havior u(t) = v(t) after transient finished. 

Noise-induced synchronization can be detected by means of direct comparison of the 
states of two chaotic systems, u(t), and, v(t), being under the influence of noise. The 
other method of diagnostics of NIS is calculating the largest Lyapounov exponent (LE) of 
dynamical system that measures the stability of the motion. Indeed, in driven chaotic system 
the largest Lyapunov exponent may become negative, that results in synchronization: both 
systems forced by the same noise "forget" their initial conditions and evolve to identical 
state [3] . If the noise influence is infinitely small the largest Lyapunov exponent is positive 
for such a system. 

In all cases of the noise-induced synchronization being known hitherto the boundary of 
the noise-induced synchronization regime is associated with the point on the parameter 
axis where the largest Lyapunov exponent of the system under study crosses the zero value 
when its sign is changed from "plus" to "minus" . In this paper we report for the first time 
that the noise-induced synchronization regime of two spatially extended uncoupled identical 
systems driven by common noise may be preceded by a new type of behavior, when the 
largest Lyapunov exponent remains zero in a finite range of parameter values. This kind of 
behavior called "incomplete noise induced synchronization" (INIS) demonstrates the features 
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of the synchronous motion of two uncoupled identical systems driven by common noise: 
although the states of the system differ, trajectories can be transformed into each other by 
an appropriate spatial shift. 

The structure of the paper is the following. In Sec. [I] we describe the system under 
study being the complex Ginzburg-Landau equations driven by common noise and the new 
type of their behavior (incomplete noise induced synchronization) observed for the certain 
set of the values of the control parameters. In Sec. [Til we consider the mechanisms being 
responsible for the INIS regime. Sec. II III deals with the different models of noise with the 
distinct probability densities. The final conclusions are given in Sec. [TV] 



I. SYSTEM UNDER STUDY AND INCOMPLETE NOISE-INDUCED SYN- 
CHRONIZATION 



The system under study is represented by a pair of uncoupled complex Ginzburg-Landau 
equations (CGLEs) driven by common noise: 

u t = u- (I- i(3)\u\ 2 u + (1 + ia)u xx + D((x,t), 
v t = v - (1 - i(3)\v\ 2 v + (1 + ia)v xx + D£(x,t), 

where u(x,t), v(x,t) are complex states of the considered systems, a and (3 = 4 are the 
control parameters, D defines the intensity of a noise term £ = Qi We have used model 

noise with the asymmetrical probability distribution of the real £i and imaginary (2 parts 
of the random variable 

p(fi 2) j2C ll2 ,ifO< Cl , 2 <l, (2) 
I 0, otherwise 

on the unit interval £1,2 G [0; 1]. The simulation of the random variables £1,2 with required 
probability distribution p(Ci,2) was carried out in the same way as it was described in 1^] 
for the exponential stagger distribution. Equation (pQ) was solved with periodic boundary 
conditions 

u(x, t) = u(x + L, t) and v(x, t) = v(x + L, t), (3) 
with all numerical calculations being performed for a fixed system length L = 40n and 



random initial conditions. To evaluate (0Q) the standard numerical scheme has been used 18| , 
the value of the grid spacing is Ax = L/1024, the time step of the scheme At = 2.0 x 10~ 4 . 
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FIG. 1: The evolution of the difference of the system states \u{x, t) — v(x, t)\ described by complex 
Ginzburg-Landau equations dU) (a) without noise and (b) with noise with the intensity D = 3. In 
the second case the difference of the states of both systems in every point of space tends to be 
zero (after transient), which means the presence of the noise-induced synchronization regime. The 
control parameter values are a = 2, (3 = 4 

If the noise intensity is equal to zero (D = 0) and initial conditions u(x, 0) and v(x, 0) are 
not identical, both systems demonstrate the complex chaotic behavior (both in time and in 
space), with the system states being different, i.e., u(x,t) ^ v(x,t) (Fig. dja). Alternatively, 
if the noise intensity D is large enough the states of both systems coincide with each other 
(Fig. [H&), that is the evidence of the noise-induced synchronization. 

To detect the presence of the noise-induced synchronization regime the averaged differ- 
ence 

t+T L 

e = -^jr J J \u(x,t) — v(x,t)\ dxdt, (4) 

T o 

between the spatio-temporal states of two CGLEs driven by common noise was calculated. 
The averaging process starts after a long-time transient with duration r = 200. 

In the NIS regime the relation e = takes place, since in this case the difference between 
the states of two identical spatially extended systems in every point of space tends to 
zero. We have also calculated the largest Lyapunov exponent A for one of the systems (JTj). 
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FIG. 2: The dependencies of (a) the averaged difference (j4]) and (6) the largest Lyapounov exponent 
of the CGLE on the noise intensity D for the different values of the control parameter a. Curves 1 
correspond to the case of a = 1, curves 2 were calculated for a = 2. The values of noise intensity 
corresponding to the onset of noise-induced synchronization are shown by arrows with labels Dj^ IS 
and D^js f° r the curves 1 and 2, respectively. The boundary of the incomplete noise- induced 
synchronization is also shown by arrow marked as Din is 

As it was mentioned above in the NIS regime the largest Lyapunov exponent A should be 
negative. 

The dependencies of the largest Lyapunov exponent \(D) and the averaged difference 
e(D) on the noise intensity D are shown in Fig. [2] for two different values of the control 
parameter a. It is easy to see that for the control parameter a = 1 (curves 1 in Fig. [2j a, b) 
the value of the noise intensity D for which the largest Lyapunov exponent A crosses zero 
value and becomes negative coincides with the point where the averaged difference (jlj) starts 
being vanishingly small. So, in this case the noise-induced synchronization boundary is 
Djsris ~ 1-5 ( see arrows in Fig. [2J a, b) and we deal with the occurrence of the noise-induced 
synchronization regime being typical and well-known. 

Alternatively, a different scenario is observed in the same system ([1]) if the control param- 
eter value a = 2 is considered (see curves 2 in Fig. El a, b). For such a choice of a-parameter 
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FIG. 3: (Color online) The dependencies of the largest Lyapunov exponent A on the noise intensity 
D calculated for the different values of the spatial grid spacing: □ — Ax = L/2 10 , At = 2.0 x 10~ 4 , 
— Ax = L/2 11 , At = 1.0 x 10" 4 , V — Ax = L/2 12 , At = 5.0 x 10" 5 . The plateau where A = 
is shown. The control parameter values are a = 2, (5 = 4 

value the largest Lyapunov exponent becomes equal to zero for the large enough intensity 
of noise Dj^js ~ 1.53 whereas the averaged difference e between the spatio-temporal states 
of two CGLEs driven by common noise exceeds zero value sufficiently (Fig. [2J a,b). With 
further increase of the noise intensity D (when D is equal to D^is ~ 2.5) the value of e be- 
comes equal to zero (see Fig. a) and the largest Lyapunov exponent starts to be negative 
which is the evidence of the presence of the noise-induced synchronization regime. 

In other words, there is the finite interval of the noise intensity values (Dinjs; D^is) for 
which the noise-induced synchronization is not observed, and the largest Lyapunov exponent 
A is equal to zero. To prove this fact we have calculated the largest Lyapunov exponent of 
the complex Ginzburg-Landau equation for different values of the spatial grid spacing. The 
range of the noise intensities corresponding to the plateau where A = is shown in Fig. [31 
One can see that the largest Lyapunov exponent calculations with the different values of the 
spatial grid step give the similar results. Based on these calculations we come to conclusion 
that the largest Lyapunov exponent is actually equal to zero in the finite range of the noise 
intensity. 

Despite the fact that the noise- induced synchronization is not observed in the region where 
A = (see Fig. [2]), this range of the noise intensities corresponds to the behavior showing 
the features of synchronous dynamics. The manifestation of synchronism can be observed if 
one of the complex media described by the Ginzburg-Landau equation starts to shift along 
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FIG. 4: The dependence of the difference s between the states of the media u(x,t) and v(x,t) 
described by the complex Ginzburg-Landau equations (UJ) on the space shift S for the control 
parameters a = 2, /3 = 4, D = 2 

the second one with the spatial shift 5. In other words, if one uses the shifted state of one 
of the system v — v{x + 8, t) in Eq. (pQ) the averaged difference e changes depending on this 
shift 5. This movement of one of the systems supposed to be very slow for the transient to 
be completed. In this case such a spatial shift 5$ may be found that both Ginzburg-Landau 
equations start to behave identically, with the largest Lyapunov exponent being equal to zero. 
Therefore, we have called this regime "incomplete noise-induced synchronization" (INIS). 

This statement is illustrated in Fig. [4] where the dependence of the difference e (J4]) on 
the space shift 5 is shown. One can see that there is value 5q of shift 5 for which the 
averaged difference e becomes equal to zero. Therefore, for this space shift Sq both systems 
demonstrate identical behavior and the noise-induced synchronization is observed. This shift 
Sq depends on the initial conditions. For the other values of the spatial shift 5 the system 
states (both in space and time) are different, but the largest Lyapunov exponent is always 
equal to zero for the considered set of the control parameter values. 



II. MECHANISMS RESULTING IN THE INCOMPLETE NOISE-INDUCED 
SYNCHRONIZATION REGIME 

Let us discuss the mechanisms resulting in the occurrence of the incomplete noise-induced 



synchronization regime. In work {4] it has been shown, that for dynamical systems with small 
number of degrees of freedom the mechanisms of the onset of noise-induced synchronization 
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and generalized synchronization are equivalent. The mechanism of the generalized synchro- 
nization occurrence can be considered with the help of the modified system approach as it 
was done in Ref. 19j for the chaotical systems with the small number of degrees of free- 
dom and in Ref. |20j for the spatially extended system. We interpret the mechanism of the 
onset of incomplete noise-induced synchronization in a similar way. Therefore, following 
Ref. {4, Ji], Q we consider the dynamics of the modified spatially extended system with the 
additional term determined by the mean value of noise. 

The deterministic modified Ginzburg-Landau equation with the additional term, deter- 
mined by the mean value (-DC) °f the noise term ( in stochastic equation (CQ) can be written 
as 

9u m , . . 1 2 / 1 . s d u m 



u 



1 - i(3)\u m \ 2 u m + (1 + ia)—^ + (DQ. (5) 



r\ , ~"ffL \^ T' / I ""ill I ""III ' \ — ' "~~V r\ o 

at ox z 
For the selected kind of noise with the probability distribution ([2]) (DQ = 2D/3. 



Equation 



e.g 
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13) is forced CGLE, widely studied and well documented in the literature (see, 
23J). It is well-known, that different types of the spatio-temporal patterns 
may be observed depending on the domain of the control parameter values. If the value 
D is large enough, the homogeneous stationary state Uq = uo(x,t) = const is observed in 
the system ([5]). In this case the largest Lyapunov exponent is negative, with the stationary 
state regime in the system (JSJ) corresponding to the noise-induced synchronization in the 
system ([1]). With decrease of the noise intensity D the stationary state Uq loses its sta- 
bility that corresponds to the boundary of the noise-induced synchronization of the initial 
Ginsburg-Landau equations (CQ) driven by noise. 

At the same time the loss of the stability of the homogeneous stationary state occurs in 
different ways depending on the control parameter values of the modified Ginzburg-Landau 
equation (JSJ). 

Indeed, the homogeneous stationary state u can be obtained numerically from equation 

M - (1 -2/3)k| 2 w + 2,0/3 = 0, (6) 



e.g., by a Newton method 24(. To analyze the stability of Eq. ([6]) we have to consider 
the linearization of the modified Ginzburg-Landau equation in the vicinity of the station- 
ary solution uq. Let u = u r + iui be a small perturbation of the homogenous stationary 
state uq = u r + iui, i.e., u m = u + u. Having linearized equation (jSJ) and assuming that 
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FIG. 5: The dependencies of the real part of the eigenvalues A on the wave number k for the 
different values of L>-parameter when the control parameter a has been fixed as (a) a = 1 and (b) 
a = 2 

u r (x, t) = u r (k) exp(At + ikx), Ui(x, t) = Ui(k) exp(At + ikx) we obtain the dispersion rela- 
tion 

1 - uf - 3u 2 - -((3u 2 r + 3{3u 2 + 
2f3uiU r — k 2 — A 2u r Ui — ak 2 ) 

[3u 2 — 2uiU r + 2f3u r iii — u 2 — 
3f3u 2 - ak 2 3u 2 + 1 - k 2 - A 
determining the stability of the homogenous stationary state Uq. The homogenous stationary 
state Mo is stable if the condition ReA(fc) < 0, VA; is satisfied. 

The evolution of Re A(/c) with the decrease of D- value for a = 1 and a = 2 is shown in 
Fig. [5j a and Fig. [51 b, respectively. One can see, that for a = 1 the homogenous stationary 
state Mo loses its stability when D « 1.5. In this case the spatial perturbation with the wave 
number k = starts growing exponentially. As a result, the stationary state uq becomes 
unstable, the spatio-temporal chaos taking place in system ([5]). With the largest Lyapunov 



0. 



(7) 
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FIG. 6: The profiles of the spatial stationary states |uo| 2 and |iifc(x)| 2 observed in ([5]) for the 
different values of the noise intensity D and a = 2. Dashed lines correspond to the unstable state 

N 2 

exponent becoming positive both in the modified (jSJ) and original (pQ) Ginzburg-Landau 
equation, the noise-induced synchronization regime in Eq. (OQ) is destroyed. 

For the value of the control parameter a = 2 the homogenous stationary state uq loses its 
stability for D « 2.5 and the spatial mode with the wave number k = ±0.5 becomes unstable 
in contrast to the case of a = 1 considered before (see Fig. [5l b). Therefore, for a = 2 the 
periodic spatial state Uk(x) = Uk(x + l) (where / is close to 2ii/k due to periodical boundary 
conditions), which is stationary in time replaces the homogenous state u in the modified 
Ginzburg-Landau equation. The example of the profiles of such stationary in time but pe- 
riodic in space states is shown in Fig. El Obviously, for such stationary states the largest 
Lyapunov exponent is equal to zero. Evidently, in the initial Ginzburg-Landau equation 
driven by noise, D((x,t), with the mean value (D() the stationary in time and periodical 
in space structure Uk{x) is perturbed by the fluctuations. Therefore, the spatio-temporal 
dynamics of Wfc(x) looks like aperiodic motion, with the largest Lyapunov exponent being 
also equal to zero. Since two identical media, u(x, t), and, v(x,t), driven by common noise 
start with different initial conditions u(x, 0) and v(x,0), the spatially periodical structures 
do not coincide with each other, i.e., Wfc(x) ^ v k (x), but there is such a shift in space 5 
depending on the initial conditions u(x, 0) and v(x, 0) where Uk{x) = Vk{x + So). Therefore, 
for Din is < D < D^is Ginzburg-Landau equations ((Tj) driven by common noise are char- 
acterized by zero largest Lyapunov exponent and their states are not identical. If the first 
of the systems is shifted along the second one with a certain shift So that, depending on 
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initial conditions, satisfies the requirement Uk{x) = Vk{x + So), the identical behavior of both 
considered systems is observed. 

Note, also, a very good agreement between the values of the noise intensity D correspond- 
ing to the loss of the stability of the homogenous stationary state u (see Fig. \5\ b) and to 
the point where the largest Lyapunov exponent becomes equal to zero (see Fig. [21 b.) 

It should be noted that the INIS phenomenon considered above is determined by the 
peculiarity of the periodic boundary conditions. Due to the usage of such a kind of the 
boundary conditions, any spatially periodic solution can be moved arbitrarily along the 
spatial coordinate, and, therefore, in this case there is additional translational degree of 
freedom. Evidently, zero Lyapunov exponent corresponds to this translational degree of 
freedom, whereas all other Lyapunov exponents are negative. Instantaneous states evolved 
from the different initial conditions under the influence of the given realization of noise can 
be transformed into each other by means of the corresponding spatial shift. Taking this 
aspect into consideration the INIS regime does not seem to be observed in the absence of 
the translational invariance, i.e., for the other kind of the boundary conditions. 

To illustrate this statement we have considered the complex dynamics of Ginsburg- 
Landau equation ([1]) with the same set of the control parameter values but with the al- 
ternative types of the boundary conditions eliminating the translational invariance in the 
system under study. We have calculated both the largest Lyapunov exponent and A and the 
averaged difference e for two types of the boundary conditions: 



du 
dx 

dv 
dx 



x=0 



2=0 



du 
dx 

dv 
dx 



= 

x=L 

(8) 



x=L 



and 

Reu(0,t) = Reu(L,t) = 1, 
lmw(0,t) = Imu(L.t) = 1. 

(9) 

Rev(0,t) = Rev(L,t) = 1, 

lmv(0,t) = lmv(L,t) = 1. 
The obtained curves A(D) and e(D) are shown in Fig. [7J 

Evidently, in both cases there is no translation invariance in the system, and, as a result, 
the largest Lyapunov exponent is close to zero, but negative, in the range of D-values 
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FIG. 7: (Color online) The dependencies of (a) the largest Lyapounov exponent of the CGLE 
and (b) the averaged difference (jl]) on the noise intensity D for control parameter values a = 2 
and (5 = 4:. Curves 1 (♦) correspond to the boundary conditions (JSj) and curves 2 (□) — to the 
boundary conditions ([9]), respectively. 

where for the periodical boundary conditions ([3]) the INIS regime is observed. The averaged 
difference e, in turn, is equal to zero in this range of the noise intensity values, that shows 
the impossibility of the INIS regime existence if there is no translation invariance in the 
system. 

III. THE INCOMPLETE NOISE-INDUCED SYNCHRONIZATION AND THE 
NOISE CHARACTERISTICS 

To illustrate that the onset of the incomplete noise-induced synchronization regime is 
caused by the mean value of noise only, we examine how the different model noise influence 
the considered media described by CGLEs (pQ). One of the typical probability density is 
the Gaussian one, therefore, it seems to be reasonable to consider the dynamics of the com- 
plex Ginzburg-Landau equations ([T]) driven by common noise with the Gaussian probability 
distribution of the real and imaginary parts of the random variable 

^ = v^ exp ("^)' (10) 
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FIG. 8: (Color online) The dependencies of the largest Lyapunov exponent A on the value a 
calculated for the noise with Gaussian distribution (I10p for the different values of the variance: 

— a = 0.1, H a = 0.2, □ — a = 0.5. The plateau where A = is observed. The control 

parameter values of complex Ginzburg-Landau equation are a = 2, (3 = 4 
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FIG. 9: (Color online) The dependencies of the largest Lyapunov exponent A on the value a 
calculated for the uniform distribution of the random variable with the different values of the 

variance: — a = 0.5, H a = 1.0, □ — a = 1.5, x — a = 2.0. The plateau where A = is 

observed. The control parameter values of complex Ginzburg-Landau equation are a = 2, (5 = 4 

Since in Sec. [J and Sec. [Ill we were not able to separate the influence of the mean value 
and variance from each other, we can do it easily for Gaussian probability density. To do 
that, the value of the noise intensity D = 1 has been fixed. In this case the mean value of 
the random process is governed by the choice of a-parameter, while the intensity of noise is 
determined by the variance a. The random variable with the required probability density 
has been generated following Ref. 24|. 

We have calculated the dependence of the largest Lyapunov exponent A on the mean 
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value a of the probability density (TTOT) for different values of the variance a. The results of 
these calculations are given in Fig. [8] where the curves A vs. a- value are shown. As it follows 
from the numerical calculations the incomplete noise-induced synchronization is observed 
in the range a £ [1.0; 1.7] (as well as in the case considered in Sec. [I]) for all values of the 
variance a. Therefore, the value of the variance a (i.e., the noise intensity) does not seem 
to play the key role in the occurrence of the INIS regime. 

Having examined the spatio-temporal behavior of two complex Ginzburg-Landau equa- 
tions ([1]) driven by common noise with the probability density (flUj) we have found for the 
range a £ [1.0; 1.7] that the states of them do not coincide with each other therefore there is 
no noise-induced synchronization as expected. Nevertheless, there is the spatial shift 5 for 
which the dynamics of system u(x,t) and v(x + S,t) is identical, which is also the evidence 
of the presence of the incomplete noise-induced synchronization regime. 

To compare the results described in Sec. [Til and Sec. Ill 1 1 with each other we have to take 
into account that the value 2D/3 has been substituted in Sec. [TTJ for the mean value of 
noise (see the explanation given below Eq. (jHJ)). Obviously, if we substitute the same value 
2D/3 for the mean value a (i.e., denoting a = 2D/3) in the case of noise with the Gaussian 
probability distribution, we obtain, that the range of the INIS-regime is D £ [1.5; 2.5] that 
agrees very well with the results of the analytical consideration given in Sec. [TTJ 

The very same results have been obtained for the uniform distribution of the random 
variables with the mean value a and variance a (Fig. [9]): there is the range of the values 
of a-parameter where the largest Lyapunov exponent is equal to zero and the incomplete 
noise-induced synchronization regime takes place. 

Based on this consideration we come to conclusion that the occurrence of the incomplete 
noise-induced synchronization regime is determined by the mean value of noise, whereas the 
variation of it practically does not play any role. 

IV. CONCLUSION 

In conclusion, we have reported for the first time a new type of noise-induced synchronous 
behavior occurring in the spatially extended systems. Such a type of incomplete noise- 
induced synchronization differs remarkably from all other types of synchronous behavior 
known so far. It may be observed in a certain range of the noise intensity values, where 
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the largest Lyapunov exponent is equal to zero and the states of two identical spatially 
extended systems driven by common noise are different, although there is an indication of 
the synchronism: if one of the systems is shifted along the second one on the certain shift the 
identical behavior of the considered systems is observed. The theoretical equations allowing 
to explain the mechanism resulting in such a type of behavior have also been given, and they 
are in perfect agreement with the numerically obtained data. The influence of a different 
model noise on the occurrence of the noise-induced synchronization regime is considered. 
Though the INIS regime has been observed here in the complex Ginzburg-Landau equations 
driven by common noise with non-zero mean value, we expect that the very same type of 
behavior can be observed in many other relevant circumstances. Since the noise influence 
may result in a pattern formation (see, e.g., 25j) we suppose that incomplete noise- induced 
synchronization can be also observed for noise with the zero mean value, with the other 
types of the spatio-temporal patterns (e.g., traveling waves) being observed. 
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